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Abstract 

A composite grid was generated in an attempt 
to Improve grid quality for a typical turbine blade 
with large camber in terms of mesh control, smooth- 
ness, and orthogonality. This composite grid con- 
sists of the C-grld (or 0-grld) In the Immediate 
vicinity of the blade and the H-grid In the 
upstream region and In the middle of the blade 
passage between the C-grlds. It provides a good 
boundary layer resolution around the leading-edge 
region for viscous calculation, has orthogonality 
at the blade surface and slope continuity at the 
C-H (or 0-H) Interface, and has flexibility In 
controlling the mesh distribution In the upstream 
region without using excessive grid points. This 
composite grid eliminates the undesirable qualities 
of a single grid when generated for a typical tur- 
bine geometry. 

A finite-volume lower-upper (LU) Implicit 
scheme can be used In solving for the turbine flows 
on the composite grid. This grid has a special 
grid node that Is connected to more than four 
neighboring nodes In two dimensions and to more 
than six nodes In three dimensions. But the 
finite-volume approach poses no problem at the 
special point because each Interior cell has only 
four neighboring cells In two dimensions and only 
six cells In three dimensions. The finite-volume 
LU Implicit scheme was demonstrated to be robust 
and efficient for both external and Internal flows 
In a broad flow regime. 

Introduction 

Turbine blades are often designed to have sub- 
stantial thickness and camber and rounded leading 
and trailing edges. For a typical turbine blade 
the H-grld does not provide a good boundary layer 
resolution around the leading edge. If very fine 
meshes are used to resolve this, the number of grid 
points In the upstream region becomes excessive. 
Thus mesh points would be wasted. A turbine flow 
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analysis on an H-grld Indicated that the standard 
H-grld needs cusps at the leading and trailing 
edges to reach a convergent solution and that It 
causes somewhat large entropy changes around the 
leading edge, where the H-grld has a large grid 
skewness and slope discontinuity. 1 These numer- 
ical errors around leading edges would be reduced 
If a C-grld were used because It provides a better 
resolution around the leading edge. A standard 
C-grid was generated for a core turbine stator 
vane. 2 It reveals large grid skewing on the 
suction side of the flow passage, which should 
Increase numerical errors. 3 * 4 For a typical tur- 
bine blade, no single grid offers satisfactory grid 
properties In the entire turbine stator or rotor 
passage. 

A two-dimensional 0/H patched grid used In a 
turbine cascade computation^ shows the slope 
discontinuity at the 0-H grid Interface and has 
extremely large-aspect-ratio meshes near the lead- 
ing edge. The numerical scheme used In that work 
Is a cell-centered scheme based on the Beam and 
Warming approximate factorization. 6 Although the 
flexibility of the scheme Is Indicated In two 
dimensions, Its likely limitations In three dimen- 
sions suggest an alternative approach. 

Two- and three-dimensional composite grids 
were generated In an attempt to Improve the grid 
quality In terms of mesh control, smoothness, and 
orthogonality. A numerical scheme that will run 
on this composite grid Is discussed. 

Composite Grid 

A stator vane ring and the vane geometry at 
mean section are shown In Figs. 1(a) and (b), 
respectively. This vane ring was used for an annu- 
lar cascade experiment reported In Ref. 2. To 
alleviate numerical errors associated with the 
H-grld skewness for the blunt leading edge, a • 
standard C-grld (Fig. 2(a)) was generated to exam- 
ine any possible advantage over the H-grld for this 
turbine vane geometry. But the grid became very 
skewed on the suction side of the passage because 
the blade was highly cambered (Fig. 2(b)). In 
addition, the control of the mesh distribution In 
the blade region was limited with the C-grld, 
because It affected the mesh distribution In the 
upstream region. Results of numerical simulations 
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using this grid are likely to be affected by the 
grid skewing. 3 

Two- and three-dimensional composite grids 
were generated as shown In Figs. 3(a) and (b), 
respectively. The blades had a constant profile 
from hub to tip and were stacked at the trailing 
edge (Fig. 3(c)). Only a selected number of grid 
lines are shown for Illustration. This composite 
grid consists of the C-grld (or 0-gr Id ) In the 
Immediate vicinity of the blades and the H-grld In 
the upstream region and In the middle of the blade 
passage between the C-grlds. The C-grld (or 
0-grld) portion can be generated by using either 
the elliptic method 7 or the algebraic method. 8 
At the C-grld and H-grld Interface (or 0-H Inter- 
face) the slope continuity was preserved so that 
no special numerical approximations are needed for 
the derivatives at the Interface. The C-grld (or 
0-grld) Is orthogonal to the blade surface and pro- 
vides a good boundary layer resolution near the 
leading edge. This composite grid has better 
smoothness and orthogonality and provides more 
flexibility In controlling the meshes than any 
single grid for a typical turbine blade with large 
camber and a rounded leading edge. 

The C-grld portion of the composite grid was 
generated by the elliptic grid generation code. 9 
Two subroutines of the code were modified. 18 The 
outer boundary subroutine was modified, to Improve 
generation of the periodic boundaries, by using the 
mean camber line for high-solidity blade rows. The 
inner boundary subroutine was modified to allow a 
more general clustering of points about the leading 
and trailing edges of highly cambered turbine 
blades. Figure 4 Illustrates the portion of the 
C-grld between n = 0 and n = n c that was 
retained for the composite grid. The grid lines 
at n = 0 and n = n ma x are lnn ^ r and outer 
boundaries, respectively. The choice of nc Is 
arbitrary and may depend on the extent of the shear 
flow region. The n 1 coordinate Is a translation 
of the n coordinate by one pitch in circumferen- 
tial direction. For the H-grld In the middle 
between the blades, the C-grld points at n = n 

c 

and n = n , on the suction side and at 
c*l 

n‘ = and n' = n^ on the pressure side are 

used In the cubic spline to preserve slope conti- 
nuity at the C-H Interface and to effect a smooth 
change in mesh size. The three-dimensional grid 
In Fig. 3(b) was constructed algebraically from the 
composite grids of the hub and shroud surfaces by 
using the cubic spline and exponential stretching. 

A two-dimensional composite grid was also gen- 
erated for the first-stage stator and rotor of the 
space shuttle main engine (SSME) fuel- turbopump tur- 
bine. Figure 5 shows the construction of the C-H 
composite grid for the rotor blade. Figure 6 shows 
the composite grid generated for the stator vane. 
One merit of the composite grid is the flexibility 
In constructing the H-grid In the upstream region. 
One can choose the number and distribution of 
meshes, and the H-grld can be extended In the 
upstream direction Independently of the C-grld part 
of the composite grid around the blade. With this 
flexibility the downstream end of the grid for the 
stator can easily be matched with the upstream end 
of the grid for the rotor In the vaneless space 
between the stator and rotor. Figure 7 Illustrates 
that the composite grid can be easily constructed 
for a stator/rotor Interaction study. 


Numerical Scheme 


Semidiscrete Finite-Volume Scheme 


The two-dimensional C-H composite grid has an 
unusual grid point that Is connected to more than 
four neighboring nodes (Fig. 8(a)). At this point 
the usual differencing techniques cannot be applied. 
If a standard finite-difference scheme were used to 
solve a flow problem on this composite grid, It 
would be difficult to treat this special point. A 
finite-volume scheme 11 to be used with this com- 
posite grid presents no problem and requires no 
special treatment because each Interior cell of 
this composite grid has only four neighboring cells 
In two dimensions and six neighboring cells In 
three dimensions (Fig. 8(b)). The finite-volume 
scheme Is described briefly here. 

The Euler equations In Integral form can be 
written as 


a_ 

at 



w dQ + 



dS = 0 


( 1 ) 


for a fixed region Q with boundary dQ. Here 
w represents the conserved quantity, F Is the 
corresponding flux term, and t Is time. 


A convenient way to ensure a steady-state 
solution Independent of the time step Is to sepa- 
rate the space and time discretization procedures. 
In the semidiscrete finite-volume scheme one begins 
by applying a semidiscretization In which only the 
spatial derivatives are approximated. To derive a 
semidiscrete model that can be used to treat com- 
plex geometric domains, the computational domain 
Is divided Into quadrilateral cells. Assuming that 
the dependent variables are known at the center of 
each cell, a system of ordinary differential equa- 
tions Is obtained by applying equation (1) sepa- 
rately to each cell. These equations have the form 


dt < s ij w u> ♦ Q u 


( 2 ) 


where S^j Is the cell area and Q-m Is the 
net flux out of the cell. This can be evaluated as 
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£ f k - ax k 9 k ) (3) 

where f^ and g^ denote values of the flux 
vectors f and a on the k*h edge, Ax^ and 
Ay^ are the Increments of x and y along the 
edge with appropriate signs, and the sum Is over 
the four sides of the cell. The flux vectors are 
evaluated, for example, by averaging the values In 
the cells on either side of the edge: 


The scheme constructed In this manner reduces to a 
central-difference scheme on a Cartesian grid and 
Is second-order accurate In space provided that the 
mesh Is smooth enough. It also has the property 
that uniform flow Is an exact solution of the dif- 
ference equations. 


L U Implicit Scheme 

The composite grid Is generated to be used for 
both Invlscld and viscous flow calculations. For 
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viscous calculations the C- or 0-meshes must be 
very fine to resolve the boundary layer. And It 
Is likely that the time step Imposed by an explicit 
stability bound will be much less than that Imposed 
by the accuracy bound of an Implicit scheme. Since 
an obvious way to accelerate convergence to a 
steady state is to Increase the size of the time 
step, an Implicit scheme is expected to have a 
faster convergence. Although the alternating 
direction implicit (ADI) scheme has been valuable 
In two-dimensional problems. Its Inherent limita- 
tions In three dimensions suggest an alternative 
approach. An LU Implicit scheme was demonstrated 
to be efficient and robust for both external and 
Internal flows in a broad flow regime. 12-14 The 
scheme was extended for three-dimensional flows on 
an H-grld. 15 This scheme will be used for the 
numerical simulation of turbine flows on the 
composite grid. 


accuracy drops to first order. The unfactored 
Implicit scheme (eq. (7)) produces a large block- ; 
banded matrix, which Is very costly to Invert and ' 
requires huge storage. An unconditionally stable 
Implicit scheme that has error terms at most of 
order (At) 2 In any number of space dimensions 
can be derived by LU factorization 

+ 0 At <D~A* * D y B + )J jl ♦ 0 At (oV 4- D y B')j 

x 4W + At R = 0, (8) 

where D x and Dy are backward difference oper- 
ators and DJ and DJ are forward difference 
operators. The reason for splitting Is to ensure 
the diagonal dominance of lower and upper factors 
as well as to make use of the built-in implicit 
dissipation. 


The conservation law form of the Euler equa- 
tions In Cartesian coordinates for two-dimensional 
flow Is 

3W d£ dG 


where W Is the vector of dependent variables and 
F and 6 are convective flux vectors: 

W = (p,pu,pv,pE) T 

£ = [pu.pu 2 + p,pvu,u(pE + p)] T (5) 

2 T 

G - C pv , puv , pv + p,v( P E + p)] 

where p, u, v, E, and p are density, velocity 
components, total energy, and pressure. The pres- 
sure Is obtained from the equation of state 

P « P(Y - 1)[e - £ (u 2 4 v 2 )] (6) 

where y Is the ratio of specific heats. These 
equations are to be solved for a steady-state 
aw/at = 0, where t denotes time. 

Let the Jacobian matrices be 

3F aG 

A = aw ’ 6 = iw 

and let the correction be 

4W » W n+1 - W n 
where n denotes the time level. 

The linearized implicit scheme for a system 
of nonlinear hyperbolic equations such as the Euler 
equations can be formulated as 

|l + 0 At (D x A + D y B)J 4W + At R = 0 (7) 

where I Is the identity matrix and R Is the 
residual 

R « D x £(W n ) 4 0 y G(w") 

Here D x and D y are central difference oper- 
ators that approximate d/ax and a/ay. 

If 0 s 1/2, the scheme remains second-order 
accurate In time; for other values of B the time 


Here A + , A' , B+, and B* are constructed so 
that the eigenvalues of " + ,, matrices are nonnega- 
tive and those of matrices are nonpositive. 

A + - l (A ♦ r A I), A" = \ (A - r A I) 

B* = j (B + r B l), B" = l (B - y) 

where 

r A > max( |\ A I), r B >max(|x B l) 

Here and represent eigenvalues of 

Jacobian matrices. Equation (8) can be inverted 
In two steps. The LU Implicit scheme needs the 
Inversion of sparse triangular matrices, which can 
be done efficiently without using large storage. 
This scheme has only two factors In three 
dimensions . 

C oncluding Remarks 

Two- and three-dimensional composite grids 
were generated to Improve grid quality for an annu- 
lar turbine cascade In terms of smoothness, reso- 
lution, and orthogonality. This composite grid, 
which has a C-type (or 0-type) grid In the Immedi- 
ate vicinity of the turbine blade, provides a good 
boundary layer resolution around the leading-edge 
region for viscous calculation, has orthogonality 
at the blade surface and slope continuity at the 
C-H (or 0-H) Interface, and controls the mesh dis- 
tribution In the upstream region without using 
excessive grid points. This composite grid elimi- 
nates the undesirable qualities of a single grid 
when generated for a typical turbine geometry. 

The C-H composite grid has flexibility In con- 
structing the H-grld In the upstream region so that 
the grids for the stator vanes and rotor blades of 
the SSME fuel-turbopump turbine can easily be 
matched In the vaneless space and can be used for 
the stator-rotor Interaction study. 

A finite-volume lower-upper (LU) implicit 
scheme Is to be used In solving the turbine flows 
on the composite grid. This grid has a special 
grid node that Is connected to more than four 
neighboring nodes in two dimensions and to more 
than six nodes In three dimensions. But the 
finite-volume approach poses no problem at the 
special point because each Interior cell has only 
four neighboring cells In two dimensions and only 
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six cells In three dimensions. The finite-volume 
LU Implicit scheme was proven to be robust and 
efficient In a broad flow regime and Is expected 
to yield accurate solutions on the improved compos- 
ite grid. 
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(A) RING SHOWING CUTOUT FOR LASER WINDOW. 


(B) GEOMETRY AT MEAN SECTION. (V/V CR = CRITICAL 
VELOCITY RATIO.) 


FIGURE 1. - CORE TURBINE STATOR VANE. 



(A) THREE DIMENSIONS. 

FIGURE 2. - STANDARD C-GRID FOR CORE TURBINE VANES. 


(B) TWO DIMENSIONS. 
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FIGURE 6. - COMPOSITE GRID FOR FIRST-STAGE STATOR OF SSME 
FUEL-TURBOPUMP TURBINE. 


ROTATION 



FIGURE 7. - COMPOSITE GRID FOR FIRST STAGE OF SSME FUEL-TURBOPUMP TURBINE. 
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